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Abstract 

We show how the Abel-Jacobi map provides all the principal properties of an ample family 
of integrable mechanical systems associated to hyperelliptic curves. We prove that derivative 
of the Abel-Jacobi map is just the Stackel matrix, which determines n-orthogonal curvilinear 
coordinate systems in a flat space. The Lax pairs, r-matrix algebras and explicit form of the 
flat coordinates are constructed. An application of the Weierstrass reduction theory allows to 
construct several flat coordinate systems on a common hyperelliptic curve and to connect among 
themselves different integrable systems on a single phase space. 
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^ ■ 1 Introduction 

r — . In the classical mechanics the arrow from the initial physical variables to the action-angle variables 

^^ ' is provided by the separation of variables and then by the Arnold construction of the action-angle 

representation pF. The motion in the opposite direction ought to allow us to construct various 
mechanical integrable systems. However, in the action-angle representation all the mechanical systems 
with fixed number of degrees of freedom are indistinguishable. To describe some particular integrable 
system one should present an explicit construction of the initial physical variables as functions on the 
action-angle variables. This mapping contains all the information about a given integrable system. 
By using variety of these mappings the different integrable models may be connected together via 
the common action-angle variables. For instance, mechanical systems may be related to nonlinear 
equations and to gauge field theory. 

As an example, investigation of the finite-gap solutions of the nonlinear problems leads to the 
introduction [|2|, |3|, |4| of analytic symplectic form il g on the Jacobian fibrations and to the definition 
of the action-angle variables on the complex space of Liouville variables. In H it is shown that 
possible obstructions to the existence of global systems of action-angle variables on symplectic vector 
bundles are a nontrivial first Chern class and the presence of monodromy at singularities. Introduction 
of the action-angle representation enables ones to consider mechanical integrable systems as systems 
associated with these variables on a torus bundle with base Ai , moduli space of complex polynomials 

2<H-1 

F(X)= HiX-Xj), (1.1) 

j=i 

and with a fiber J(C), the g-dimensional complex Jacobian of auxiliary curve C defined by the Abel- 
Jacobi map U p], 0] . The fact that action-angle variables could be used for quantization of classical 
systems leads to introduction of semiclassical geometric phases. This approach results, for instance, 
in a quantum conditions on the moduli of ^-dimensional Jacobi varieties m . 

By using this Abel-Jacobi map hi and the Jacobi problem of inversion, the so-called root variables 
{Pji 9j}?=i (§7 h] on an associated Riemann surface C may be constructed instead of the action-angle 
variables. In these root variables on the level of integrals of motion the action is represented as a 



sum of items depending on one coordinate only, i.e. these variables are separated variables. The 
corresponding Ricmann surface C depends on parameters (moduli), parameterizing the moduli space 
M. of C ||, M. In terms of mechanical integrable systems the curve C is interpreted as a time- 
independent spectral curve, integrals of motion are some specific coordinates on the moduli space M 
and Jacobian J(C) is a common level of the involutive integrals of the system H. 

In what follows, we have to describe appropriate mechanical systems together with their phase 
space in initial physical coordinates {pj, Xj}™ =l . In particular, separated coordinate systems ought to 
be orthogonal curvilinear coordinate systems on the flat Riemannian manifold mU O. In this case, 
these separated coordinate systems are associated to some solutions to the Lame equation jl^, [l4], [l5| . 
Recently, the solutions to this equation have been obtained in an explicit form with the help of the 
"dressing procedure" |0|, the Baker- Akhiezer function 14 Hq] and the Lie algebraic construction 
jfq] within framework of the inverse problem method. 

The main objective this paper is to illustrate how fixed mapping from the action-angle variables 
|Q. [| Q to separated variables completely defines all the principal properties of mechanical systems. 
We shall consider the uniform Stackel models associated to the Abel- Jacobi map U on the hyperelliptic 
curve C and the well-known elliptic, parabolic and spherical curvilinear coordinate systems on IR . 
Also we discuss relations of these mechanical systems with other integrable models associated to the 
same algebraic curve. 

2 The Stackel systems 

One of the oldest problem of the hamiltonian mechanics is to find the quadratures for the integrable 
hamiltonian systems. The simplest models integrable in quadratures are the Liouville systems and 
the Stackel systems JlTj ( the Liouville systems are a particular case of the Stackel systems). 

Before proceeding father it is useful to recall the classical work of Stackel [17] . The system asso- 
ciated with the name of Stackel ||17| is a holonomic system on the phase space IR , their hamiltonian 
is 

n 

H = Y fSj {q 1 ,...,q n )(p} + U j ) . (2.1) 



i=i 



,2)1 



Here {Pj,9j}? = i are canonical variables in IR with the standard symplectic structure and with the 
following Poisson brackets 

n 

Q„ = ^2 d Pj A d 1j > {Pjilk} = 5jk ■ (2-2) 

There is an even stronger version of the Stackel theorem. 



Theorem 1 For a hamiltonian system with hamiltonian H of the form (2.1) the following assertions 
are equivalent: 

1) The associated Hamilton- Jacobi equation is separable. 

2) There exists a nondegenerate n x n Stackel matrix S, whose elements 5kj depend only on qj 

det S ^ , — -^ = , for j^m 

oq m 

and such that (2-3) 

n 

^2%kj(<lj)gj(qi, ■ ■ ■ , q n ) = hi ■ 

3=1 

3) There exist n functionally independent integrals of motion which are quadratic in momenta. 

Let C = [cik] denotes inverse matrix to S such that Cj\ = g,-. Then the Stackel theorem JlT], O 
asserts that there are n first integrals of motion, namely 

n 

h=Y J ^{p 2 3 +U ] ) , h=H. (2.4) 



The common level surface of these integrals 

M a — Iz e R n : Ife(z) = a*,, fc = 1, . . . ,n\ 
is diffeomorphic to the n-dimensional real torus and one immediately gets 



PI 



OS 



dq 3 



-7T- I = 53 "fe- 5 ^ (ft ) ~ ^ (ft ) . 



(2.5) 



fc=i 



where S(qi . . . , q n ) is an action function Q]. If this real torus is a part of complex algebraic torus 
then the corresponding mechanical system is called an algebraic completely integrable system |18 



The Stackel theorem allows to reduce the solut ion of the equations of motion to a problem in 
algebraic geometry. We can regard each expression (2.5) as being defined on the Riemann surface 



C, : y) = F S (A) , F j (A) = £ a k 5 kj (A) - U 3 (A) , 



(2.6) 



fc=i 



which depends on the values a k of integrals of motion. All the pairs of variables (j>j,qj) lie on 
these Riemann surfaces (2.6). Considered together, they determine an n-dimensional Lagrangian 
submanifold in R 



C {n) : Ci(pi,qi) x C 2 {p2,q2) X ••• xC„(p„,g„). 

The associated Hamilton- Jacobi equation 



dS dS dS 

"a7 + "■ (*' a — ' ■ • ■ ' a — j91j • 



0. 



TJj ds_os_ = E 

d qj d qj 



(2.7) 



(2.8) 



on the local manifold V n with diagonal metric g n = gAqi, . . . , q n ) analytic in the local coordinates 
{qj} has the following additive solution 



S(qi ...,q n )=^2 Sj(<Zj) , s 3 {q 3 ) = / yFjiq-j) dq 3 



(2.9) 



with the functions Fj(X) defined in (2.6). Coordinates qj{t,a\, . . . ,a n ) are determined from the 
equations 



$ij(\)d\ 



piA»(»,?o) \/YJk=i a kS\j{\) - Uj{\) 



= ft = t , 



(2.10) 



$hj(X)dX 



" pj{pj,ij) 
j=1 J~fo(po,qo) V52k=i akSkjW - Uj(X) 



Pk 



k = 2, ... ,n, 



where points Jj(pj,qj) and 7o(poj9o) be on the curve Cj (|2.6|). Notice, that bounded motion in this 



case will not be periodic in general but only conditionally periodic 
turning points determined by the conditions that functions Fj(X) (Ej6 
motion Wjk are equal to 

'■ A * S kj (X)dX 



w k] 



a« VFJW 



Hi 0. If A° and Xj are the 
2j| ) vanish, the periods of the 

(2.11) 



Thus, Stackel ITT]] showed that the orthogonal coordinates {qj}l—i permit separation in the Hamilton- 
Jacobi equation (2.E) if the metric 



ds 2 = 2J gjj(?i) • • • , Qn) (dq J ) 2 
i=i 



(<7i,...,g„) = gJqi,...,q n ) 



(2.12) 



is in the Stackel form 

Sjj(Qu--->Qn) = H|(<7i,...,9„) = -^jr , (2.13) 

where S J means the cofactor of Sji in matrix S ( |2.3| ). Here g,-„ is a diagonal metric and Hj are 
called the Lame coefficients. The modern approach to construction of the Lame coefficients see in 

@ @ 0- 

Henceforth, we shall restric t ou r attention to the uniform Stacke l systems , where all the potentials 



Uj(qj) = U(qj) and curves Cj (|2.6| ) are equal. Variables {sk,Wk} ( |2.9| , |2.1lD on a single curve C are 
the action-angle variables for the uniform Stackel systems. To constru ct the metric g,(gi, . . . , q n ) and 
the potentials U(qj) in an explicit form we shall iden tify periods Wk ( 2.11 ) with periods of the Abel 



differentials on a common hyperelliptic curve C (|2.5| ) alon g the elements of a homology basis 

In this case definition of the separated variables {<£,} ( |2.10 ) leads to the Jacobi inversion problem. In 



{ Qj } (ETo 

-ixSfJ) 



the next Section we prove that the Stackel matrix S ( |2.3| ) is completely defined by the derivative of 
the Abel- Jacobi map U on C at generic point (so-called Brill-Noether matrix) . 

3 Uniform Stackel systems and algebraic curves 

To begin with let us briefly recall some necessary facts about the action-angle variables on the Jacobian 
J{C) [0, o, H. The main ingredient of this construction is a universal configuration space, which is 
the moduli space |1(| of all algebraic curves with fixed jets of local coordinates at a fixed number of 
punctures. This concept is closely related to the notion of the Baker- Akhiezer function on admissible 
curves |14J and to the theory of algebraic completely integrable systems |18[ . 

Let us consider a genus g Riemann surface C with N ordered punctures Pj and with two special 
Abelian integrals y and A with poles of order at most I = (lj)^ =1 and m = (mj)^L 1 at the punctures. 
The universal configuration space M g (l, m) can then be defined as a moduli space of C under certain 
constraints on the set of algebraic geometrical data H, [1. In this case the space A4 g (l,m) is a 
complex manifold with only orbifold singularities. To introduce the local coordinates on M g (l,m) 
we cut apart the Riemann surface C along a homology basis Aj, Bj j — l,...,g with canonical 
intersection matrix 

A i oA j =B i oB j = 0, A i oB j = 5 ij . (3.1) 

By selecting cuts from Pi to other Pj for each 2 < j > N one gets a well-defined branch of the 
Abelian integrals y and A. Among the complete set of local coordinates on Ai g (l,m) the following 
moduli are distinguished 

ydX, j = l,...,g. (3.2) 



The universal configuration space A4 g (l,m) is a base space for a hierarchy of fibrations C^ k \l,m) of 
particular interest to us. These are the fibrations whose fiber above each point of A4 g (l, m) is the fc-th 
symmetric power S k (C) of C. This fiber C^ k \l,m) is the set of all effective divisors D = 71 + • • - + 7fe 
(the 7j's may not be mutually distinct) of degfc of C, i.e. C^ k '(l,m) can be identified with the set of 
all unordered fc-tuples {71, . . . , 7^}, where 7j's are arbitrary elements of C. 

Let T> be the open set in Ad g (l, m), where the zero divisors of dy and dX do not intersect. Fixing 
all the local coordinates on Ai g (l,m) except Sj ( [3.2] ) one can determine a smooth p-dimensional 
foliation of V, independent of the choice we made to define the coordinates themselves Q. Hereafter, 
by abuse of notation, one leaf of this foliation, is denoted just by M. and C^ means the above 
fibrations restricted to M.. 

Let dS = ydX be a meromorphic 1-form on C with the special Abelian integrals y and A, which 
have fixed expansions near the punctures Pj |4[] . It means that we have imposed the certain constraints 
on the algebraic geometrical data (according to |L4J we used admissible data). These constraints 
ensure the existence of global system of action-angle variables and the presence of the corresponding 
symplectic form |g| . The fact that we impose some constraints provides us with additional properties 
of dS. Namely, generating 1-form dS possesses the property 

ddS dydX . 

^- = ^7 = dW >> 3 = 1,-- -,9, 0.3) 



where Sj are action coordinates ( |3.2| ) on the moduli space M and differentials of the angle variables 
dwj form some basis of holomorphic differentials (normally, even if differential is holomorphic, its 
moduli-derivative is not). Moreover, form dS give rise to differentials spanning a whole space Tt\(C) 
of holomorphic differentials. Hence, for any generic divisor D = 71 + • • ■ 7 g on C the standard 2-form 
on C^ 

(9 \ g a 

is a desired holomorphic symplectic form fl g on C*- 9 '. The set of variables {s.,-, Wj} 9 , 1 are the complete 
set of action-angle variables on J(C). These action-angle variables {$j,Wj} have been obtained by 
generalizing the definition of actions introduced for integrable systems on tori in the form of periods 
of holomorphic differentials dwj along the elements of a homology basis in EL p| . 

Now we turn to the uniform Stackel systems. The corresponding algebraic curve ( p.q ) is a 
hyperelliptic curve given by an equation of the form 

2g+l 

and puncture P is the point at infinity A = 00. Recall that the moduli Xj of C are integrals of motion 



(2.6). Solution to the inverse Jacobi problem and associated Abel-Jacobi map on C relate a set of the 



action-angle variables and the separated variables. 

Variables of separation qj(t) give solution to the inverse Jacobi problem (2.1C). The associated 
Abel-Jacobi map U : Div(C) — > J(C) is restricted to Lagrangian submanifold C^ 

U: C (k) ^J{C). (3.6) 

Note that whenever we discuss the Abel-Jacobi map, we shall tacitly assume that a base point 70 



( 2.10) on C has already been fixed in an appropriate position. 

Suppose that point D = 71 + • • • + 7fc, k < g belongs to C^ k \ The differential of the Abel-Jacobi 



map (3.6) at the point I? is a linear mapping from the tangent space Tn(C^ 9 ') of C^ at the point D 



into the tangent space T U ^(J(C)) of J(C) at the point U{D) 

W D : T D {C^)^T U{D) {J{C)). 

Now suppose that D is a generic divisor, and Zj is a local coordinate on C near the point jj. Then 
(zi, . . . , Zk) yields a local coordinate system near the point D in C^ k >. Let dwk (k = l,...,g) is a 
basis for a space Tii(C) of holomorphic differentials on C, and near jj 

dw k = (f>kj(zj)dZj , 

where ipkj(zj) is holomorphic. ft follows that the Abel-Jacobi map U can be expressed near D as 

U(z x , ...,z k )= 1^2 / <t>ij(zj)dzj, ...,^2 / <t>gj{zj)dzj 

W=l ^70 j—i J 10 

Hence 

01l(7l) ••• 4>gl(ll) 

w D =\ ■■ •■. ; I • (3.7) 



(t>lk{lk) ■■■ 4>gk(lk) 



is the so-called Brill-Noether matrix. 



Theorem 2 Transpose Brill-Noether matrix Up on the genus g > n hyperelliptic curve C, which is 
the derivative of the Abel-Jacobi map IA at generic divisor D , deg-D = n, is equal to the Stackel matrix 
S for the uniform Stackel system on C with metric 

det S 
gjj{qi,...,q n ) = —j -■ 



At generic point D , degD = g matrix S = Up is regular matrix satisfying t he S tackel theorem. 

At g > n we have to consider restriction of the Abel-Jacobi map (3.6) onto C^ n K In this 



case symplectic form fl n on C^ n ' is an appropriate projection of Q g (|3.4[) and C'™) be a Lagrangian 



submanifold in the phase space C . The separated variables {Pj, 9j}?=i are constructed from the 
first 2n action-angle variables (^9|) only and the action differential dS = Y^^iPj^Qj §i ye rise to an 
n-dimensional chart of the whole space TLi(C). The corresponding n X n Stackel matrix is the left 
upper nxn block of the general matrix S 
n = g. 

As an example, let us consider some basis for Tii(C), for instance 



Up and, therefore, unless otherwise indicated, we assume 



dwi 



x>- 



y(A) 



dX, 



.7 = 1, 



,9- 



(3.8) 



By choosing this basis we fix a basis of action-angle variables (3.2-3.4). To solve the Jacobi inversion 
problem (2.1C) one gets variables of separation 



Pj=y(lj), %=A(7i)» i = i,- 



,9 



(3.9) 



for a generic point D = 7i + • • • + 7 S on C, which coincides with divisor of simple poles of the 
corresponding Baker- Akhiezer function p4| . In the real case (when p j an d qj are real), the separated 
variables qj (3.9) (so-called root variables) vary along cycles Aj (3.1) over basic cuts on C and, 
therefore, our problem is defined on g-dimensional real torus. The holomorphic symplectic form Q g 



on C^ coincides with standard ones (2.2) and a fiber C^ be a complex Lagrangian submanifold of 
the phase space C 9 



(^ZD 

C (n) s gn(Q . 



(C(A) x C(fi) x ■ • • x C(y)) /a n 



n< 3: 



(3.10) 



where a n is the permutation group on n letters. 

Recall, that derivative Up bears a great resemblance to the usual Gauss mapping. The map 

U* D induces a canonical mapping from C into the (g — l)-dimensional projective space C — > P 
On the other hand, the canonical mapping is defined the derivative of the Abel-Jacobi map. For a 
hyperelliptic curve C of genus g > 2, the canonical map C — > P 9 is the composition of the double 



covering map C — ► P , sending (y, A) to A, with the Veronese map P — ► I 
the polynomial ring of degree g — 1. With respect to the basis of Hi(C) (3 
C has an extremely simple expression 

By using this map we introduce the g x g matrix 

/ X^ 1 



i)9-l 



given by a basis for 
the canonical map of 



,A 9 " 



,A,11 



,.s-i 



,9-i 



S(A,M> 



\ 



A 
V 1 



1 



v 
1 / 



(3-11) 



determined on a Lagrangian submanifold (3.10). For a generic point D = 71 + 
Stackel matrix is equal to 

s kj ( qj )= x°- k 



j g in (3.£) the 



S(qi,q2,.-.,q g ) = S(A,/x, 



W=qi,p,=q2,—,u=q g 



A=(jj 



(3.12) 



Recal l, that the diagonal metric g ? - 3 is completely determined by the corresponding Stackel matrix 
( 2.13 ). Nevert heless , we introduce another eq uiva lent definition of the metric. Substituting the 
Stackel matrix ( 3.12 ) in the algebraic equation (2.2) one gets 



^2$kj(qj)gjj(qu(to, ■■■,q g ) = Ski 

3=1 



3=1 



x k ~ l 

•» e(A) 



(3.13) 



1 r a^ 1 

2?ri J c e(A) 



where by definition 



Sjj(Qi,Q2,---,q n ) = Res| ; 



e(A)' 



Here we introduced function e(A), which has zeroes at the points qj giving solution of the inverse 
Jacobi problem. 

In general, function e(A) with g zeroes, which are solution of inverse Jacobi problem, is expressed 
in the Riemann theta-function 



e(A)=0(W( 7 i,...,7 s )-/3-#) 



W(7i 



t7g )=U( 7 i) + ---+U( Jg ). 



Here K is a vector of the Riemann constants and j3 — (/3i, . . . , j3 g ) £ 
principal properties of the function e(A) (3.14) are considered in ||. 



is a fixed vector 



(3.14) 



The 



Proposition 1 Function e(A) on C with g zeroes (p j , qj) giving solution to the Jacobi inversion 
problem is completely defined the metric gjj(qi, q2, ■ ■ ■ , Qn) l ift- IS ) for a uniform Stackel system. 



We prove this proposition in the polynomial ring only. In this case 



e(X) = l[(X-q k ), 



and 



,(<?!, 92,- 



k=l 



.«») = Res \x= qj e '(A) = g -pr 



To prove ( 3.13| ) for this metric, it suffices to consider the following integral 



ce(A) 



77- <P "TTT = V Res| 



='i.i 



3=1 



e(A) 



= - Resl 



e(A) 



= S k 



k,g-l 



(3.15) 



(3.16) 



(3.17) 



where C encloses all qj. 

Function e(A) is defined on the universal configuration space, i.e. it is independent on the moduli 
Xj of C (integrals of motion) and on a choice of the basis of holomorphic differentials in Tii(C). For 
instance, in the polynomial ring let us consider a set of the equivalent Stackel matrices with the 
following entries (3.7) 



SkjWL = <f>kj(X)\, 



''1.1 



cj) kj (X) = XV- k +a^XV- k - 1 + ... + a^ k , 



(3.18) 



where polynomials <fikj form various basises for the polynomial ring of degree g— 1. Substituting (3.18) 
in (3.13) and (3.17) one obtains at once universal solution e(A) (3.15). Below we shall see that the 
hamiltonian H (2.1) with the diagonal metric g ? „ (3.16) is closely related to the distinguished puncture 
P at infinity A = oo on the hyperelliptic curve C (3.5). The different Stackel matrices (3.18) correspond 
to the distinct sets of the integrals of motion in the involution for a single hamiltonian H. The 
completeness and functional independence of these integrals directly follows from the completeness 



and independence of the basis elements (3.18) for a polynomial ring 



Finally, we look at other fibrations C^ n> at n ^ g. At g > n, to construct the metric gjAq 



on C^ n \ we expand the initial curve C (3.5) by 



2g+l 



2n+l 



= n ( a - a >) = u *9+iw + n ( a - a >) < n ^ a ■ 



,q n ) 



(3.19) 



Here U2 g +i(X) is an at most 2g + 1 order polynomial, which is regarded as a potential in (2.6). The 
n x n Stackel matrix and the corresponding function e(A) may be associated to the auxiliary genus 
n curve 



2n+l 



C: 



~2 

V 



]\{X-Xi). 



(3.20) 



Function e(A) are independent on the moduli of C 



3.19) and, therefore, uniform potential U2 g +i in 



(3.19) has an arbitrary form and decomposition ( [3. IS ) determines the highest power of the polynomial 
U(X) only. 

At n > g the above holomorphic symplectic form Q, g on the leaves M. is degenerate. However, a 
non-degenerate form on &- n > may be obtained by restricting C^ to the larger leaves M. of the foliation 
P|. The leaves M. correspond to the level sets of all the local coordinates except to holomorphic sj 



(3.2) and to some additional (n — g) coordinates associated to meromorphic differentials duij in (3.3- 
3.4). In fact, to construct the action-angle variables we have to add several meromorphic differentials 
to holomorphic angle variables. Thus, at n > g the symplectic 2-form J7„ on C'"' is meromorphic B. 
As an example, at n = g + 1, we can add one local coordinate in the neighborhood of puncture 
P at infinity B. This additional coordinate occurs in the Stackel matrix and in the metric in the 
following way 



S( ff+1 >(A, M ,...,^) = goS (9) (A,/v..^) 



;(A) = q a JJ(A - qj) g 00 = Res| A=c 

3=1 



(3.21) 



A 9 " 



5(A) 



At n > g these systems with meromorphic form Q n po ssess several reductions of the additional 
meromorphic coordinates, for instance go = const in ( 3.21 ) |20(| . 

Above formulas are well adjusted for generalization. If the curve C (B^q) is substituted by 



C 



y 2 = F(\) 



PKA) _ ng 1 (A-A J ) 
QmW nr=i(A-4) 



m < 2g + 1 , 



where {5k} is a set of m arbitrary constant, one gets 



•Sfej(A)|; 



e(A) 



n^i(A-g,-) 

Qm(X) 



lA= ^ Qm(X) 

Note, that the algebraic equation ( 3.13| ) is covariant with respect to the transformations 

S^i?- 1 (A)S, e(A)^i?- 1 (A)e(A), 



(3.22) 



(3.23) 



that leads to the general form of the metric 



gjj(qi,Q2,---,q n ) = Rcs| 



9m(A)fl(A) 

m =1 (A-%) 



(3.24) 



associated to the curve C. We shall use this freedom to consider the standard curvilinear coordinate 
systems @ H H- 

So, the hyperelliptic genus g curve C may be associated to a family of the uniform Stackel systems 
on the phase space C by using the Abel-Jacobi map U, its differen tial IA*^ and their restrictions 
on C". Diagonal metric g.-.-(<7i, 92, ■ • • ,Qn) ( 2.13 ) in the hamiltonian (2.1) is completely defined by 
number of degrees of freedom n and potential U(X) is at most 2g + 1 order arbitrary polynomial. 

On the other hand, one fixed metric gjj(qi, ■ ■ ■ , q n ) may be associated to an infinite set of the 
hyperelliptic curves C. The corresponding hamiltonian systems differ from each other by the power 
and by the form of polynomial potentials U(X) ( p. 19 ). Among these systems we must to distinguish 
systems on C^ n > at n > g ( p.10 ) for which the number of degrees of freedom n is more than genus g 
of the associated curve C. In this case the corresponding symplectic 2-form on C is meromorphic Q. 
In the next section, we shall identify these systems with the degenerate or superintegrable systems 
pl[ . Recall, that for degenerate system the number of independent integrals of motion is more than 
number degrees of freedom. 



4 The Lax representations 

Let us recall that the key idea, which has started the modern age in the study of classical intcgrablc 
systems, is to bring them into the Lax form. All the properties of the uniform Stackel systems may 
be recovered from the properties of the Abel map. Nevertheless, now we want to obtain the Lax 



representations for all the uniform Stackel systems associated to the hypcrclliptic curve C(y, A) (3.5). 
We consider construction of the Lax representation as a necessary intermediate step to study quantum 
counterparts of the Stackel systems. 

In the simplest case the Lax matrices L{\) or L(y) are defined as the matrix valued functions 



on bare spectral curves F\, X £ F\ (1.1) or F v , y G F y , while the full spectral curve C(y, A) is given 
by the Lax eigenvalue equations 

C: det (L(X) - y) = , det (L{y) - A) = . (4.1) 

As a result, C arises as a ramified covering over the bare spectral curve F\ or F y |22[ . 

Till now a delicate questions is how to construct the Lax matrices L(X) or L(y) for a given 
integrable system. The one integrable system may be associated to the different curves and one curve 
C may be associated to the different mechanical integrable system on a common phase space. As an 
example, the ro-particles Toda lattice can be equivalently formulated in terms of two different Lax 
representations p3| associated to the single hyperelliptic curve C. 

Here we consider equation for a general algebro-geometric symplectic structure associated to the 
spectral curve C of the given Lax representation L 

Q n = - V Res Pa v — (4.2) 



proposed in g. Here f2 n is the restriction of the algebro-geometrical symplectic form ( p.4|) on C 
generated by two differentials dy and dX having poles at punctures P a . Functions ip and ip + are the 
Baker- Akhiezer function on C and it's dual function. If we fix some 2-form fi n and the Baker- Akhiezer 
functions ip, ip + on a given curve C, then one can attempt to recover the associated Lax matrix L. 

For the some particular Stackel systems the 2x2 Lax matrices [E(J |24j and the corresponding 
vector Baker- Akhiezer function ip associated to natural vector fields on the Jacobian of any hyperel- 
liptic curve are known. On the other hand, we know the general scalar Baker-Akhiezer function if) 
on C defined by its analitical properties on C, which corresponds to geodesic systems with diagonal 
metric @. 

Note, here we have the vector Baker-Akhiezer function ip, which is the eigenfunction of the matrix 
L associated to the curve C, and scalar Baker-Akhiezer -0, which is completely defined by analitical 
properties on the same curve C. 



For the uniform Stackel systems let us identify the preassigned symplectic structure Q (2.2) with 



the symplectic structure (3.4) defined on a hyperelliptic algebraic curve C. Next we try to recover 
Lam matrix for a geodesic motion under the following additional assumptions: 

1) L(X) is a generic 2x2 matrix associated to a spectral hyperelliptic curve C of genus g — 
[(n-l)/2]. 

2) The associated vector Baker-Akhiezer {unction ip has a constant normalization a J& 

< a, i\) >= aiV>i + a 2 ip2 = 1, a = (0, 1) . 



3) The first component of 'tp i n (J^.i ) is proportional to the unique Baker-Akhiezer function ip on 
C with fixed analytical properties jj^|/. 

At the first assumption n is a number of integrals of motion, which are moduli of C ( n = 1g + 1) 



and, therefore, form fl n in ( |4.2| ) is a restriction of meromorphic symplectic form fl g (3.4) to the 
minimal n-dimensional leaf M. |^ for integrable systems on C . The second ass ump tion allows us 
to reduce vector Baker-Akhiezer function to scalar one. In this case, solution of ([4.2| ) is completely 
defined by the function ip on C only. At first we present this particular solution in term of the function 
e(A) associated to the Abel map IA. Introduce function e(A) and its time derivative 

e(A)= rc^-ft) m <„ ; e x (X) = {H,e(X)}, (4.3) 



where {6 k} is a set of m arbitrary constant and H be a hamiltonian of the geodesic motion 

1 






,Qn)p 



% jj {qi,...,q n ) = Res| A=?j . -^ 



(4.4) 



Thus, in the Lax equation for a geodesic motion 

L X (X) = {H,L} = [L,A] , 
matrices L and A are given by 



L(X) = 
A{\) = 



-ex/2 e 



(A) 



/i e 
/ -ft 



(A). 



(4.5) 



1 




The hamiltonian H (4.4) is equal to a highest residue at the distinguished Weierstrass point on C at 
infinity A = oo 

H = - Res| A=DO A"-" l detL(A) . 



where the full spectral curve C is equal to 

C: y 2 = F(X)=dctL(X), 

F(X) = -h 2 (X)-e(X)f(X) = 



(4.6) 



(4.7) 



2 4 

By definition zeroes of e(A) be separation variables and conjugated variables pj are given by 

j(qi,...,q n )Pj 



Pj = MA)L 



-ij 



h{\) = -ex/2 






A-9i 



(4.8) 



In accordance with p5| pairs of separation variables (qj,Pj) lie on the spectral curve C 



y 2 ( lj )=P 2 =h 2 {X)\ x= 



\=qj 



■F(X = Qj ) 



F(X)\ 



1i 



As usual, rational function F(X) admits some different representations 

n fc =i(A-4) ^(a-4) fc= ^ +1 



(4.9) 



Here {/,}" =1 and {Jk}k=i are two sets of independent integrals of motion in the involution. The 
first set of integrals {Ij} 7 I = i in ( |4.9D corresponds to the Stackel matrix ( 3 . ll| ) . The set of equivalent 
Stackel matrices fl3.18j ) relate to another decompositions of a numerator of F(X) ( ff.9| ). The second 
set of integrals { Jk}k=i m (4-9) is associated to an expansion near punctures {5k, oo} on C. 
The spectral curve C (|4.l|) is a time-independent curve and, therefore, 



{H,F(X)} = 0, 



d 3 x e(X) 



= 0. 



(4.10) 



Thus, in fact |2q] , we consider the polynomial solutions e(A) = J1(A — 9j(£)) to the equation ( 4.10 ) 
and describe the hamiltonian dy nam ics of their zeroes qj (t) (re call, th at d x means derivative by time) . 
Substituting fun ction e(A) (13) and hamiltonian H ( |4.4| ) into ( 4.10| ) one gets the equations in 



the metric g 7 -.. = H 2 ( 2.12 ). If we introduce so-called rotation coefficients 



A- 



d.M 3 



i^j, 



(4.11) 
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these equations may be reduced to the following equations |13J 

dkPij = PikPk 3 , i^ j ^k, 

difaj + 8jf3j l + E PmiPmj =0, i ^ j, 



(4.12) 



where the notation i =/= j =/= k means that indices i,j, k are distinct. 

Of course, these equations may be obtained without any Lax representation by using definition 



(2.13) of the metric, properties of the Abel-Jacobi map and preassigned asymptotic behavior of e(A) 
at the distinguished point A = oo. 



The equations (4.12) are equivalent to the vanishing conditions of all a'priory non-trivial compo- 
nents of the curvature tensor JRl O, FL5). Therefore, using (4.12) we conclude that local Riemannian 
submanifold (V„, g| v ) (2.i) of the Riemannian manifold (C",g) is a flat manifold whose metric is 



diagonal with respect to the coordinates {qj}. Imposing some additional restrictions on the space of 
solutions to (fO)J2(l, one could get the Bourlet type equations |t5| related to another Riemannian 
manifolds of constant curvature. 



To construct more general solutions to ( |4.2| ) associated to hypcrelliptic curve C of higher genus 
we begin with calculation of the Poisson bracket rela tions for the initial Lax matrix L(X). It allows 
us to identify the space of solutions to equation (L2) with the loop algebra £(sl(2)) in fundamental 
representation |23| and then to use the representation theory of the underlying algebra sl(2) J27). 

Theorem 3 The Poisson bracket relations for the matrix L(X) ( ETa J are closed into the following 
r-matrix algebra at m < n only 



{L(A),2(/i)} = [ri 2 (A,/x),i(A)] - [rai (A, /*),!(/*)] 



(4.13) 



Here the standard notations are introduced: 

L(X) =L{X)®I : 

n 



fl2 (A, /i) 



A — jJL 



L(p) = I®L{n), 



r 2 i(A,/u) = ILri 2 (fJ-, A)II. 



(4.14) 



and II is the permutation operator of auxiliary spaces P3| . 

The Poisson bra cket relations for the Lax matrix L(X) (4.5) are preassigned by the initial sym- 
plectic structure (3.4). It is necessary to calculate two brackets only 



{ e (A),e(/i)} = 0, 



(4.15) 



and 



{h(X),e(p)} 



*WE 



g ■ (qi,..., q n )Pj E[?=i ( A - Qj ) 



3=1 



X-qj 



nr=i(A-^) 



'(A)e(M)E 



~i (A - qj)(fi - qj) 



(4.16) 



e(A)e(/z) 



E 



A — [i *— ' \ A — qj fi — qj J X — fi 
where we used a standard decomposition of rational function 



[e(p) - e(A)] 



e- X (A)-E 



i=i 



( L, 



Resume" 1 (A). 
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Another Poisson bracket s m ay be directly derived from these brackets and by definition of the entries 
of the Lax matrix L(X) (|4.5|) via derivative of the single function e(A) 



{h(\),h(n)} 
{/(A),e(/x)} 

{/(A), h^)} 
{/(A), /(/*)} 



0, 
9 x {ft(A),e(/x)} 



A — jjb 



MX) - h(n)] 



(4.17) 



1 



-\dl{h{\) 1 e{ii)} = 
2 A — fi 

- l -dl{h{X) 1 e{^)} = Q 1 



[/(A)-/(M)], 



To derive the first bracket we have to combine second and first derivatives of the brackets (4.15) and 
(4.16), respectively. At the last bracket one substitutes the equation of motion ( 4.10J ). 

If, contrary to our geometric conventions, the order of polynomial Q m (A) is more then order of 
polynomial Pi in ( 3.22| ), i.e. if m > n in the metric (4.5), then rational function e(A) admits another 
representation 



-i(A) 



3JJ 



S*"* 



+ £(\,qi,...,q n ), 



where remainder £(A) is a certain polynomial. Substituting this function e(A) into (4.15-4.17) one 
gets 

9£(A, qi,...,q n ) 



dx 



= o. 



This constraint to remainder £(\,qi, ■ ■ ■ ,q n ) directly follows from the symmetry of the last Poisson 
bracket in ( j4.17p . 

The r-matrix algebra ( 4. 13| ) is so-called linear case of the r-matrix algebras corresponds to 
intcgrable systems, which are modelled on coadjoint orbits of Lie algebra s/(2). The r-matrix in 
(4.14) is a standard rational r-matrix on C(sl(2)) |28| . The general form of the function e(A) (3.14) 
leads to the elliptic and trigonometric r-matrices ]28|, |2Sfl 



Thus, for a geodesic motion (4.4) the Lax representation ( |4.5| ) with arbitrary poles {5k}™=i (4.3) 
may be regarded as a generic point at the loop algebra C(sl(2)) in fundamental representation after 
an appropriate completion j28[ . Sin ce, to construct the Lax representation for a potential motion 
with the fixed metric g._- (gi, . . . , q n ) ( LS ) we can use the outer automorphism of the space of infinite- 
dimensional representations of sl(2) proposed in fl27f . 

Applying this automorphism of the underlying algebra si (2) directly to the Lax representation 
L(A) ( |4.5| ) on C(sl(2)) we obtain a family of the new Lax pairs 



L'(A)=i(A)- ( r_.[0(A) e - 1 (A)] 



JV 





1 



A'(A )= ^-a_.[0(A )e - 2 (A ) ] JV= ( u ; (A) J) 



(4.18) 



Here </>(A) is a function on spectral parameter and \z]n means restriction of z onto the ad^-invariant 
Poisson subspace of the initial r-bracket |27], |2^, Q ■ For the rational r-matrix ( 4.14 ) we can use the 
linear combinations of the following Taylor projections 



Z N 



+ OG 



E ^ 



,fc= — CX3 



JV 



N 

E 

fe=0 



z k \ 



(4.19) 



or the Laurent projections [gT], p9J] . 

The mappings (4.18) from the representation of the loop algebra £(sl(2)) to representations of 
the universal enveloping algebra U(C) play the role of a dressing procedure allowing to construct the 
Lax matrices L' N (A) for an infinite set of new integrable systems starting from the single known Lax 
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matrix L(A) associated to one integrable model. This mapping preserves the metric g,((?i, . . . , q n ) in 
the hamiltonian (2.1), but changes the potential U(qj) and associated curve C. 

New Lax matrix L'(X) (4-18) obeys the linear r-bracket (4.13), where constant r^-matrices 
substituted by r\ -matrices depending on dynamical variables |27[ [29 1. 



r\i{\\i) -* r' 12 =r 12 



([<f>(X)e~ 2 (X)} N -[^)e- 2 ( f ,)} N 
(A-M) 



(4.20) 



We have to distinguish systems on C^ at n > g ( [3.1C| ) for which the number of degrees of freedom 
n is more than genus g of the associated curve C. According to J4| the corresponding symplectic form 
is meromorphic. In this case the action differential dS = ydX give rise to a whole space TLi(C) and, in 
addition, several meromorphic differentials on C. We can identify these systems with the degenerate 
or superintegrable systems [gl| . 

Theorem 4 The complete set of noncommutative integrals of motion for the degenerate uniform 
Stdckel systems with meromorphic symplectic form Q g is determined by the generalized spectral surface 



C{y,X,n) 



det (yl + ILL' (A) ® !/(;«)) = . 



Here we used the outer product of the 2x2 Lax matrices L'(A) with L'(fi) and II means 4x4 
permutation matrix in C x C . Equation of motion for the matrix L(A, /j) = IIL'(A) £g>L/(/x) is equal 
to 



-L(A,/z) = L(X, f i)A(X^)-UA(X, f i)U- 1 L(X, f i), 
dt 



(4.21) 



A(X,fx) 



A(X) ® I + 1 ® A(p) 



where matrix A(X) is a second Lax matrix and / is a unit matrix. 

It is easy to derive from (4.18), that n > g iff n > N, where N is a highest power in the Taylor 
projection (4.19). In this case the corresponding r-matrix (4.20) preserves the simple pole at the 
puncture P at A = oo and the associated second Lax matrix A 1 remains a constant in spectral sense 
dA{X) 



dX 



under the mapping (4.1 



Thus, for the degenerate systems A(X, /i) = IL4(A, /i)n _1 and equation (4.21) takes the standard 
Lax form and it proves the theorem. 

As usual, spectral curve C ( |4.l| ) of L'(A) is a generating function of the involutive family of 
integrals of motion. Substituting functions 0(A) = A™Q~ 1 (A)L T /v(A) into L'(A) ( |4.18| ) one gets their 
spectral curve in the form 



C 



y 2 = F'(X) = dctL'(A) = U N (X) + 



m=i(A-<y' 



where {/'■} are integrals of motion. It is a time-independent curve and, therefore, 



dF'(X) 
dt 



= 0, 



-dl + u N (X)d x + -ujv,x(A) 



e(A) =0. 



(4.22) 



Of course, this equation may be obtained directly in framework of symplectic geometry [^lj). Let us 
briefly explain an origin of this equation in the theory of nonlinear equation, that allows us to relate 
scalar Baker- Akhiezer function ip and function e(A). 

The same algebro-geometrical symplectic form fl g (B-4 ) on hyperelliptic curve C ( |3.5| ) leads 
directly to a hamiltonian structure for soli ton equations |2|, |4j]. As an example, we consider the KdV 
equation associated to hyperelliptic curve ( |3.5| ) with one puncture P (N — 1) at infinity A = oo and 
at I = 1, m = 2 H. Let us select one leaf of foliation corresponded to dX with all zero periods 



dX = 



13 



for an arbitrary cycle C. In this case, the Abelian integral X(P) is a single-valued function, with only 
a pole of second order at P (m = 2). For finite-gap solutions of the KdV equations, moduli Sj ( |3.2[ ) 
are canonically conjugated with respect to the Gardner-Faddeev-Zakharov symplectic structure to 
angle variables Wj (see |32| and references within) . Thus, the uniform Stackel systems have a common 
set of the action-angle variables with solutions of the KdV equations. 



Starting with this set of variables we consider general algebro-geometric equation (4.2) for non- 
linear systems. Solution of the equation (4.2) in a ring of second order differential operators with the 
standard Baker- Akhiezer function ip on C d3.5| ) is well known M, H, EL M . The associated Shrodingcr 
operator has the form 

C(X)=--^+u(x,t,X), (4.23) 

where A is a parameter. In some simple cases, such as the KdV equation, this parameter A appears 
as an eigenvalue and one ultimately equates the potential u with a solution of the nonlinear equation 
itself. Let us look for a solution .4(A) of the Lax system in the ring of differential operators 



£(A)^ = o, 



(4.24) 



f dC(X) 

\ at 



[C,A] V = 



of the form 



.„> ,.,3 1 de(X) 
Substituting the given form of A into the Lax system, one gets 



(4.25) 



du _ 

at ~ 



1 -dl + u{X)d x + \u x {X) 



e(A). 



(4.26) 



Equation ( 4.26 ) is called the generating equation. For a different choices of the form of e(A) and 
u(X), this procedure leads to different hierarchies of integrable equations, as an example to the KdV, 
nonlinear Shrodinger and sine -Gor don hierarchies pL |9| or to t he D ym hierarchy [p3| . If we consider 
the solutions of the equation ( 4.26 ) in the form of polynomial ( 3.15 ), then the roots qj of e(A) define 
the root variables and as a result finite-gap solutions of the problem of geodesic (see || |^, |33| and 
references within). 



Substitution of the special form of second operator A(X) ( 4.25 ) into the Lax system ( 4.24 ) allows 
us to eliminate the Baker- Akhiezer function ip and to construct 2x2 Lax matrix in e(A). In fact, we 
replace the Baker-Akhiezer function ip on C to the mutually disjoint function e(A) on C, which has a 
transparent mechanical interpretation (3.15). Recall, that function e(A) is defined as function with 
zeroes, which give solution to the Jacobi inversion problem H on the hyperelliptic curve C. 



5 The flat coordinates 



According to p4[ at n = g the orthogonal curvilinear coordinates {Pj,Qj}^ = i form a generic divisor 
of the simple poles of the Baker-Akhiezer function ip, which is defined by their analytical properties 
on C. The evaluation of -0 at a set of punctures on C determines the flat coordinates {pj, Xj} 3 - =1 for 



the diagonal metric ( 2.1 2j ) . It turns out that up to constant factors the Lame coefficients H., are equal 
to the leading terms of the expansion of the same function ip at the punctures on C E3 . 

Next we reach the same conclusions by using the function e(A) and the corresponding Lax rep- 
resentation L(X) on C. As usual, we reduce the study of algebraic geometrical data to the analysis 
of the associated geodesic motion. The crucial observation is that the equations of motion in coor- 
dinates {pj, Xj} r f =l on the Riemannian manifolds of constant curvature have a Newton form and the 
corresponding hamiltonian has a natural form 



= £j(xi, ■ ■ ■ ,Xn) , H = ^dijpipj + V(xi, . . . ,x„) . 



eC. 



(5.1) 
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(5.2) 



where £j(xi, . . . ,x„) and potential V(x\, . . . , x„) are functions on coordinates only. Let us introduce 
new function B(X) 

S 2 (A)=e(A) = H- 2 (A), 

which is "inverse" to the Lame coefficients Hj ( |2.f3| ). One immediately gets 

>(A)- 



F{\) = B A B X 



F'(X) = B 3 B XX + B i 



B 4 



(5.3) 



A ! 



These equations have the form of Newton's equations for the function B 

B xx = F(X)B- 3 , 



(5.4) 



B xx = F'(\)B- 3 -B 
To expand function B(X) at the Lourent set 

N 



4>W 



B 4 



N 



J=0 



it is easy to prove that coefficients Xj obey the Newton equation of motion (5.4) (see ( 4.10J) and 
references within [ft6| ). Here we reinterpret the coefficients of the bare curves F(X) and F'(X) in (5.4) 
not as functions on the phase space, but rather as integration constants. In variables Xj mapping 
( 4.18 ) affects only on the potential (x-dependent) part of the integrals of motio n Ik- The kinetic 
(momentum dependent) part of Ik remains unchanged. So, the dressing mapping (4. IS) allows us to 
get over from a free motion on C to a potential motion on C . 

As an example, from (5.4) we get some well known orthogonal curvilinear coordinates on M. 
(see gjj |l| [[§): 

elliptic coordinates m — n in ( |3. 23] ) 

= l + 2^ T"^ = B ( A > x i' ■ ■ ■ ' x ») 



e(\,qi,...,q n ) 



nLi(A-4) 



k=l 



X- Sk 



Si < x 1 < S 2 ■ ■ ■ < S n < x n 



e(X,qi,...,q n ) 



parabolic coordinates m = n— 1 in (3.23) 

n- =1 (A 



on 



= x - x n + y^ 



k=i 



A - Sk 



= B (A,xi,. ..,x n ) 



xi < (5i < x 2 • • • < <5 n -i < x„ 



e(A,g , 



spherical coordinates m = n + 1 see (3.21) 

<?oirj =1 (A-^ - +l 



> 9nJ = 



ni:;(A 



E 

fc=i 



Y 2 



A 



= B (A,xi,.. .,x„ + i) 



Curvilinear coordinates {qj} are zeroes of function e(A) and flat coordinates {x^} are residues of 
e(A) = B 2 (X) at the punctures, in accordance with the Baker- Akhiezer function approach |1J, Pq| . 

All the separable orthogonal curvilinear coordinate systems in M may be obtained from these 
coordinate systems [[19], ^2] EfJ. According to [Q, all the possible separable in these coordinates 
potentials, which are polynomials or rational functions of the cartesian coordinates Xj, belong to 
the set of the uniform Stackel systems. Thus, we can claim that every such mechanical system is 
embedded into a proposed scheme. 
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5.1 Quasi-point canonical transformations. 

In conclusion, we discuss another parameterizations of the function e(A). Of course, function e(A) 
admits various representations in different variables and we can use this freedom, as an example to 
solve equations of motion |26| . The considered above parametrization describe the point canonical 
transformation only. Here we discuss an application of the Weierstrass reduction theory to construct 
another cartesian coordinates on C. 

It is obvious, that the Lax representation 

L(X) = [L(X),A(X)] 

is covariant with respect to the transformation of the first Lax matrix 

L(A)-f^(A,A 1 ,...,A fc )L(A) 

with an arbitrary function </>(A, Ai, . . . , Xk) on time-independent moduli {A^} of C and on spectral 



parameter A. However, this transformation drastically changes the Poisson bracket relations ( 4.14 ) 
and parameterization of L(A) in the flat coordinates {pj,Xj}. Hence, in addition to considered above 
flat coordinates {pj , Xj }™ =1 , the same function e(X,qi,...,q n ) may be associated to another set of flat 
coordinates. Now we show that to introduce these new variables {pj,Xj} we can use various covering 
of the initial curve C, as an example, covering listed in [ j35| ]. 

Let us assume that the initial torus J(C) = T 2g may be decomposed in a direct product of several 
tori 

T 2g = T 2 9l x . . . x T 2g k ^ ^ g . = g _ (55) 

The corresponding Riemann matrix has a block form B = B± x B% • • ■ x B/~, where Bj are the gj x gj 
Riemann matrices and the corresponding Baker- Akhiezer function on C is factorized. In this case we 
can consider curve C as a if-sheeted covering of tori T 2gj . Such covers are known to exist for any 
K > 1 and for arbitrary tori 139] 

First of all, we can introduce the separated variables {qj} associated to a whole torus T 2g . For 



dynamics on J(C) = T 2g the corresponding Lax representations L(X) (4.5) are 2x2 matrices 



Secondly, we can introduce another set of separated variables {qj} associated to each torus T 2g - 



in (5J5). For dynamics splitting on several tori T 2gj the Lax representations have a block form 

L(A)= ■, (A), (5.6) 

where Lj{X) are the 2x2 matrices defined by functions ej(A) on the each torus T 2gj pjj. Two 
sets of variables {qj} and {qj} are related by canonical transformation induced by the covering, that 



allows us to get 2x2 Lax matrix instead of matrix (5J3). It means that we have two isomorphic 
integrable systems with different Lax representations and the corresponding canonical transformation 
is a quasi-point transformation |p6| . 

To illustrate this construction we take as an example several systems at n = 2. Starting with an 
hyperelliptic curve C of genus g = n = 2 we define variables (pi,qi) and (^2,92) on the Lagrangian 
submanifold C^ 2 ' (|3.9| ). The Jacobi inversion problem is the problem of finding these variables from 
the equations ( 2.10 ) with the Stackel matrix S given by ( 3.1 2| ) . This problem is solved by using the 



Kleinian p-functions, which are second logarithmic derivatives of the Kleinian c-function 

dlna((3 u (3 2 ) 



dfrdfr 



P22 = qi + 92 , P12 = -qiq2 



(for detail see [B5L p7|). The function e(A) (3.15) on C with zeroes at the points qi, qi is equal to 



e(A) = A 2 - p 22 A - P12 = (A - qi )(X - q 2 ) = A 2 + 2A Xl + (2x 2 + x 2 ) , (5.7) 

1G 



Here we used the free dom ( 3.23 ) and cartesian coordinates {xj} or {x'j} are derived from the "inverse" 
Lame function 6(A) ( |5.2| ). Applying the outer additive automorphism of sl(2), we can construct the 
Lax matrices L'(X) for an infinite set of integrable mechanical systems with the following hamiltonians 



H = piP2+V N (xi,x 2 ), 



(5.9) 



H 



p'l+p'l + V^{x[,x 2 



Among them, we distinguish the Hcnon-Hciles systems at N = 3 and the systems with quartic 
potential at N = 4. For these systems the genus of associated curve C is equal to the number of 
degrees of freedom g = n = 2. 

Function e(A) ( 3.13J) is independent on the moduli of C and, therefore, the above construction of 
the integrable systems ( |5.9| ) readily gets over on the reducible curve C. To construct this reducible 
curve, let us take two tori T 2 2 



= £(i-0(i 



fc±0, 



(5.10) 



with a Jacobi moduli 



4 = - 



(l-a)(l-/3) 

Making the rational order two (K = 2) change of variables 



w± 



V(l-«)(l-/3) 



A=F \foiJi 



(x-a) 2 {\-py 



:V- 



(A-a)(A-/3)' 



A. 



one gets hyperelliptic curve 



C 



y 2 = A(A - 1)(A - o)(A - /3)(A - a/3) , 



(5.11) 



(5.12) 



which gives a two-sheeted covering of two tori T 2 2 ( 5.10J ). It is a well-known example of the reduction 
of hyperelliptic integrals to elliptic ones by using the rational change of variables proposed by Legendre 
and generalized by Jacobi p5| . 

The complex torus T 2 is isomorphic to the curve of genus g = 1 given by equation w 2 — /(£). In 
the above, we have presented the covering for the two odd curves (5.10) at deg(/) = 2g + 1 = 3. All 
computations concerning the even curves at deg(/) = 2g + 2 — 4 give similar covering 139], so we do 
not present these formulae. The odd and even curves at g — 1 are associated to the integrable cases 
of the Hcnon-Hailes system and system with quartic potential, respectively. 

Next we can introduce two pairs of variables (pi, q{) and (p 2 , §2) being on the tori T 2 2 . Functions 
e i,2(A) on T 2 2 are equal to 

e 1 (A) = A-q 1 , e 2 (A)=A-q 2 . (5.13) 

Variables {pj,qj} are separated cartesian coordinates for the integrable systems on T 2 x T 2 with the 
hamiltonians 

#3,4 = pI +p\ + v 3A (qi) + ^3,4(92) , (5.14) 

which is a sum of two one-dimensional hamiltonians on T 2 2 - The corresponding 4x4 Lax represen- 
tation has a block form (5.6), whose blocks are determined by the functions ei^A) (5.13). 

The covering ( 5.11 ) induces canonical transformation of variables {pj , qj } to {pj , qj } |37| . These 
pairs of variables lie on the different curves T 2 2 and C, respectively. The common moduli a and 
(3 of these curves are integrals of motion. On the orbit O (a = const, /3 = const) this canonical 
transformation (5.11) becomes a point transformation. It is so-called quasi-point transformation p6j. 

By using change of variables induced by covering (5.11) one can construct the 2x2 Lax matrix 

for the evolution (5.14) splitting on two tori. In variables {qj} matrix L(X) is determined by the 

function 

(A-a)(A-/3), 



e(A) 



-e(A), e(X) = (X-q 1 )(X-q 2 ) 



(5.15) 



(l-a)(l-/3) 
In fact, we add two additional zeroes a and /3 into the function e(A) (3.13|) on the reducible curve C 



(5.12) and, therefore, change parameterization of the Lax matrices in flat coordinates {pj,Xj}. 
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^ 2 9j 



In general, to introduce new flat coordinate s, we can take any tori T x £ of arbitrary genus 
<?i.2 > 1 and consider two-dimensional evolution (5.14) splitting on these curves with an arbitrary 
one-dimensional potentials V^+ife)- The standard change of variables 



<H 



xi ±x 2 



e(A) = A 2 - xi + 



(5.16) 



preserves the natural form of the hamiltonians (5.14) for arbitrary potentials V2 gj +i(qj). The equa- 
tions of motion remain the Newton equations in these variables {Jjj , x, } 

In the considered above example (5.1C) both in depe ndent hyperelliptic integrals are reduced to 
elliptic ones by using a common substit ution £ — > A ( 5.11 ). It relates to existence of the second order 
automorphism of a hyperelliptic curve (5.12) p5|: 



(A, y) 



a/3 

T' 



\-> 



VaZfi 



(5.17) 



It allows us to introduce another parameterization of the function e(A) in cartesian coordinates, which 
preserves the natural form of the hamiltonian. Namely, in addition to ( p. 16 ), we can use the following 
canonical transformation of variables {Pj,qj} to the cartesian coordinates {pj,Xj} 



e(A) = A 2 



Q+ + Q- , (Q+-Q-) 



xi 



4xi 



(5.18) 



Here functions Q±(pj, Xj) are the classical counterparts of the supercharges in two-dimensional SUSY 
p6[ with the following properties 

{H,Q±} = ±f(p j ,x j )Q ± , {tf,Q+Q-} = 0. 



At g = 2 (N = 3 or N — 4 in ( 5.14)) these functions Q± and / on variables {pj,Qj} or {pj,Xj} are 
listed in [p6| . Moduli a and (3 in (5.17) are integrals of motion, therefore, automorphism r induces a 
second quasi-point transformation ass ociat ed to t orus T 2 x Tf ( |5.10 ). 

Two quasi-point transformations (5.11) and (5.18) for the physical variables {xj}, {xj} and {xj} 
bind together all the integrable cases of the Henon-Hailes system at N = 3 and three integrable 
cases of the system with quartic potential at N = 4. Of course, these systems have a common set of 
action-angle variables. Moreover, the same variables are associated to the Kowalewski top j|, which 
is a supersymmetric quantum model as well. 

Thus, several supersymmetric models are related to evolution splitting on the tori, when the 
number of degrees of freedom is equal to the genus g = n of the associated covering curve C = T\ X T2. 
It would be interesting to get a geometrical interpretation of these supersymmetric objects arising 
from finite-dimensional SUSY quantum mechanics. 



6 Conclusion 

It is known, that curves y + y^ 1 = F(X) together with the 1-forms 



V 



dS & = log A^, 

y 



are implied by integrable models of the Toda chain family (standard and relativistic models). The 
corresponding Lax representations are defined on the Poisson-Lie groups with quadratic r-matrix 
algebra. The corresponding mapping from the action-angle variables to separated variables has been 
proposed in p2[ . 

On the other hand we can consider the umbilic solutions of the KdV equation P, p3| . These 
systems are defined on a generalized Jacobi variety of the symmetric product of n logarithmic Rie- 
mannian surfaces in place of the Liouville tori. Nevertheless, it is possible to introduce variables that 
linearize the corresponding hamiltonian flows. These systems may be interpreted as counterparts of 
the discrete-time Stackel systems. 
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Both these sets of models are associated to the change of parameterization of hyperelliptic curve 
from "plane" parameterization to "annulus" ones (A — > log A). The crucially interesting lift to the 
interpretation of A as a coordinate on elliptic curve. 

For all these integrable models it would be interesting to estimate the possibility of application of 
the usual Stackel approach. On this way we should consider mapping between action-angle variables 
and separated variables, and should study the differential of this map. In the presented paper there 
are the Jacobi inversion problem and differential of the Abel-Jacobi map. 

This research has been partially supported by RFBR grant 96-0100537. 
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